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Abstract 

We derive in detail four important results on integrals of Bessel functions from which three 
combinatorial identities are extracted. We present the probabilistic interpretation of these iden- 
tities in terms of different types of random walks, including asymmetric ones. This work extends 
the results of a previous paper concerning the Darling-Siegert interpretation of similar formulas 
emerging in the analysis of random flights. 
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1. Introduction 

The analysis of random nights (see Orsinghcr and De Gregorio (2007)) performed by a particle 
changing direction at Poisson paced times with uniformly distributed orientations of the steps lets 

* Corresponding author. 

Email addresses: valentina.cammarota9uniromal.it (V. Cammarota), aime.lachal9insa-lyon.fr (A. Lachal), 
enzo.orsingher9uniromal.it (E. Orsingher). 



emerge a number of intriguing relationships involving different types of Bessel integrals. These are 
the source of identities with curious combinations of Gamma functions which suggest interesting 
probabilistic interpretations. 

The instrument to give a probabilistic background to these combinatorial formulas are random 
walks on the set Z (with upward and downward steps) or on N (with non-decreasing sample paths). 
We arrive at different forms of Darling-Siegert expressions where the distribution of first-passage 
times and transition functions are suitably combined. We work on generating functions (instead 
of the reflection principle as in De Gregorio and Orsingher (2007)) and obtain directly many 
distributions, including the first passage time through an arbitrary level a when steps take values 
±1 with probability p and q. The case where some conditional first hitting time of random walks 
is related to combinatorial identities is also examined below. 

In the analysis of a planar random motion with uniform orientation of the steps, the semigroup 
property of the Bessel functions (1.5) yields as a byproduct the identity (1.9) which suggests an 
interpretation based on a non-decreasing random walk with steps taking values or 1. 

Our aim is to make the paper as self-contained as possible and for this reason all the Bessel 
integrals are proved in detail as well as the distributions used in the underlying probabilistic 
interpretation. 

We recall the usual representations of Bessel functions (Gradshteyn and Ryzhik (1981), formulas 
8.402 and 8.411.1): for \i > 0, the function is defined by the following equivalent equalities 

w-f l3 &(-r <") 

n=0 v ' 

M x ) = ir- I e ixsin0 -^ d6. (1.2) 

2?r J -IT 

We denote by N the set of the positive integer numbers including zero and by Z that of the integer 
numbers. 

Our analysis starts with the following integrals of Bessel functions. 
Proposition 1 For any [M,v,a> 0, we have that 
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Actually, the above formulas hold for larger domains of the parameters /x and v (see Gradshtcyn 
and Ryzhik (1981)), but we shall only deal with positive /i and v. 

Let us emphasize that (1.3), (1.4), (1.5) and (1.6) are particular cases of the general family of 
Bessel integrals J Q a x a (a — J^(x) J u {a — x) dx and simple formulas exist only in very few cases. 

From Proposition 1 we deduce the following formulas for the Gamma function. 
Theorem 2 For any fi, v > and r e N, we have that 
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Formulas (1.7) and (1.8) come respectively from (1.3) and (1.4) while (1.9) comes from (1.5) or, 
equivalently, (1.6). 

Formulas (1.7) and (1.8) lead to the following Darling-Siegert-type relationships (see Darling and 
Siegert (1953) and also De Gregorio and Orsingher (2007) for the discrete-time version) involving 
the distribution of the classical random walk {S n ) neN on Z with steps taking values +1 or —1. 
Let us introduce T a = min{n > 1 : S n = a} for any a e Z, with the usual convention that 
min(0) = +oo. The random variable T a is the first hitting (discrete) time of level a by the random 
walk. 

Corollary 3 For any fi,v,r £ N \ {0}, we have that 
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Formulas (1.7) and (1.8), and then (1.10) and (1.11), have already been stated in De Gregorio and 
Orsingher (2007). However, therein, they were related to the symmetric random walk instead. It 
is remarkable that all random walks with steps ±1 lead to the universal identities (1.7) and (1.8) 
whatever the probability of jumps. 

Similarly, let {S' n ) neN be the random walk with steps taking values or 1 and let T' a = min{n > 
1:5^ = a}. Formula (1.9) leads to both following Darling-Siegert-type relationships. 
Corollary 4 For any (i,i/,r£ N, we have that 
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It is noteworthy that (1.12) and (1-13) come from the only universal identity (1.9). This fact hinges 
on a feature of the non-decreasing random walk (<S^) neN and will be explained in Subsection 3.4 
and Remark 12. For the same reason, (1.5) and (1.6) are equivalent. 

Of course, formulas (1.10), (1-11), (1-12) and (1.13) can be derived by means of independent 
arguments by using the Markov property of random walks as explained in Subsection 3.1. We have 
chosen to prove them by means of Theorem 2 (or equivalently of Proposition 1) because of the 
importance of Bessel funtions in many fields of Mathematics and especially in Probability Theory. 
There is a very huge amount of formulas concerning the special functions of Mathematical Physics, 
and we find it interesting to provide for certain ones, when possible, a probabilistic interpretation. 

Our paper is organized as follows. In Section 2, we recall some well-known formulas of Bessel 
integrals. Some proofs of them relying on hypergeometrical functions can be found in Watson 
(1952), chapters 12 and 13. Nevertheless, elementary proofs are not so accessible in the literature 
and then, because of their interest, we provide them in order to make the paper self-contained and 
easily readable. In Section 3, we derive the Darling-Siegert-type interpretation of formulas (1.7), 
(1.8) and (1.9). We point out that the Darling-Siegert-type relationship appearing in Darling and 
Siegert (1953) deals with continuous Markov processes whereas our work concerns discrete Markov 
processes. 

The main tool of our analysis is the use of generating functions for the probabilities appearing 
in Corollary 3. 

2. Mathematical Background 

Let us recall well-known formulas which will be used throughout the paper: 
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2.1. Proof of Proposition 1 



The first step for deriving formulas (1.3) and (1.4) is the evaluation of the Laplace transforms of 
the Bessel function J v which is presented in the next lemma (see Gradshtcyn and Ryzhik (1981) 
formulas 6.611.1 p. 707 and 6.621.1 p. 711). 
Lemma 5 For a,B,v> 0, we have that 
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For integer values of the index v we are able to present a proof of these formulas. 



PROOF. 

• Proof of (2.4). We consider the integral representation (1.2) of the Bessel function which we 
rewrite as J v (x) = ± jZ n e ^ xsin0 - iv6 &6 = e^" 088 ^ 9 d9 and we get 

\ e- ax JJ(3x)dx = / e iv6 d6 / e -(«-tf«»«)* dx = / dQ 
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For any integer number v 1 the function f(z) = a _j /3cosz is periodic with period 2it and its poles 
are ± ^ +iargsinh^ + 2ktr, k e Z. We take a contour r# consisting of Ti = {(x,y) : x e 
[0,2n],y = 0}, T 2 = {(x,y) : x = 2ir, y e [0,i?]}, T 3 = {(x,y) : x e [0,2tt],|/ = i?} and T 4 = 



{(x, y) : x = 0, y e [0, i?]} with i? such that encloses the unique pole 8 = ^ + i argsinh j with 
residue p = r-f . = e ^fjy^i +^ ") gy Cauchy's residue theorem we obtain (2.4). 
• Proof of (2.5). It is straightforward to derive (2.5) from (2.4) as follows: 
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Proof of (2.6). In order to prove (2.6), we use the equality 
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which can be obtained by a direct calculation. Set U v = J Q °° e~ rxx x v J v (Bx) dx. In light of (2.8), 
after an integration by parts, we obtain that 
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A further integration by parts and formulas (2.8) and (2.9) together with xe~ ax = —^{{^ + 
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where we made use of the relationship J v (j3x) — 2 ^~ 1 - ) J„_i(/3x) — J v - 2 (Px) which can be proved 
by observing that ^zii _ j^^) = £~ Q - l) (f ) 2 " + ^ 2 = 

J v (px). From (2.10), (2.1) and Uq = , 1 , we extract the recurrence relationship 
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• Proof of (2.7). From (2.9), we derive J °° er ax x v J^^x) dx = | f °° e-"^ J„(/3x) dx which, in 
light of (2.6), yields (2.7). 



We now employ the Laplace transforms (2.4), (2.5), (2.6) and (2.7) to prove Proposition 1 which 
plays an important role in our paper. 

PROOF OF PROPOSITION 1. 

• Proof of (1.3) and (1.4). We can prove formulas (1.3) and (1.4) by taking the Laplace transforms 
of both members or directly as shown in De Gregorio and Orsingher (2007). Let us evaluate the 
Laplace transforms of the left-hand sides of (1.3) and (1.4). Remarking that we are dealing with 
convolution products, we easily get, by Lemma 5, 
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As a result, formulas (1.3) and (1.4) immediately emerge by inverting the Laplace transforms. 
• Proof of (1.5) and (1.6). We can prove formulas (1.5) and (1.6) in the same way. Indeed, e.g., for 
(1.5), 
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and thus, by inverting the Laplace transform, we conclude that (1.5) holds. The proof of (1.6) is 
quite similar. 



2.2. Proof of Theorem 2 



We now proceed to prove Theorem 2. 
• Proof of (1.7). We start from (1.3), using the series representation (1.1) and solving the integral 
J Q a x 2k+tl ^ 1 (a — x) 2l+,J dx, we can write that 
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By comparing the coefficients of the entire series (2.11) and (2.12), we immediately obtain formula 

(1.7). 

• Proof of (1.8). Analogously by (1.4), we write that 
r a J^(x)J v (a- x) 
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By comparing the coefficients of the entire series (2.13) and (2.14), we obtain (1.8). 
• Proof of (1.9). Similarly, from (1.5) and (2.1), we have 
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which proves, by translating [i and by 1/2, formula (1.9). 

3. Darling-Siegert-type relationship for random walks 

Formulas (1.7), (1.8) and (1.9) can be interpreted by means of random walks as we explain below. 

3.1. Background on random walks 

In order to facilitate the reading and because of its importance in many hclds, we recall some 
well-known results on random walks on Z (see, e.g., Feller (1968) Chapter 3) together with their 
proofs. Fix p € (0, 1), q = 1 — p and let (5' n ) n6N be a random walk on Z starting at the origin 
(i.e. 5*o = 0) with parameter p. For any j € Z, we introduce the generating function G(£,j) of the 
probabilities P{S n — j}, n e N: 



G(£,j) = X>{S„=j}f\ 



n=0 

Remarking that the random variable {S n +n)/2 has a binomial distribution with parameters (n,p), 
we have, for \j\ < n such that n + j is even, that 

nSn=j}= (n^)^^. (3-1) 

and this probability vanishes when n + j is odd. Notice the following equality which will be used 
further: 



nSn = -j} = [ - p ) P{S„ - j}. (3.2) 



Lemma 6 For |£| < 1 and j e Z, the function G(£, j) is given by 
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PROOF. We can directly calculate G(£, 0) and G(£, 1): since P{S n = 0} vanishes for odd values 
of n, we have, by using (2.3), that 

G(£,0) = £P{S„ = 0}r = E ( t JW 2 " = (! - « 2 )- 1/2 - 
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Analogously, since P{S n = 1} vanishes for even values of n, we have, by using again (2.3), that 



a< { , i) - £ ns„ - ir . f (* + ,y - 1 K 2 :) - ^ ( 2 ;) 



n=l n— 



n+1 „n^2n+l 



y/l - Apq£, 2 2q£ v ' 2^ v /] _ 4^ 

Now, from the elementary recurrence relationship G(£,j) — p£ G(£,j — 1) + G(£,j + 1), which 
comes from the obvious identity F{S n +i = j} = pP{S n = j — 1} + qP{S n = j + 1}, we are able 
to prove (3.3) by induction. First, for j > 2, by supposing that (3.3) holds for G(£,j — 2) and 
G(£,j — 1), we have that 

G(U) = -^G(U-1)-Zg(M-2) 

J- 2 
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Therefore, (3.3) holds for any positive integers j. Finally, we can easily check, by using (3.2), that 
(3.3) holds also for negative values of j: in fact, we can see that G(£, —j) = (^)- J G(£,j). 

Set now T a = min{n > 1 : S n = a} for any integer a with the usual convention min(0) = +oo. For 
a^O, the random variable T a is the first hitting time of level a for the random walk (S n ) n efi an d 
To is the fist return time to level for (S , Il )„ e N- 

Proposition 7 (Darling-Siegert formulas) If a and b are both positive or both negative integer 
numbers, we have that 

n — b n — b 

P{S n = a + 6} = ]T P{T Q = k} F{S n - k = 6}, P{T a+b = n} = ^ P{T Q = k} P{T b = n-k}. 
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PROOF. Let a > 1 and b > and let us consider a path connecting level at time to level a + b 
at time n. This path must necessarily pass by level a between times and n, that is, T a < n. So, 
from the strong Markov property of random walks, the translation invariance, and the fact that 
Sr a = a (in words, discrete continuity), we have 
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The same discussion holds also in the case a < — 1 and b < 0. 



Lemma 8 The generating function of T a is given by 
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PROOF. From Proposition 7, we have, for a > 1, 6 > or a < -1, b < 0, that = a + &} = 

XX=o P{^o = ^} P^n-ife = o} which involves a convolution product. So, we derive 
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G(Z,a + b) = J2nSn = a + b}C= ^P{T a = nK" J>{S n = 6}f 
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= E(£ T -l {To<+oo} )G(£,6) 



and then, for 6 = 0, E (£ T » l {Ta<+oo} ) = §§^y, which yields, thanks to (3.3), result (3.4). 



Corollary 9 The distribution of T a is given as follows: for 1 < \a\ < n such that n + a is even, 
HT a = n} = M (Z^P^^ = H P{S n = a} (3.5) 

and if n + a is odd, this probability vanishes. Moreover, the probability of eventually hitting level a 
is given by 
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PROOF. From (3.1), (3.3) and (3.4), we note that for a > 1, 
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By integrating this equality with respect to £, we obtain 



E(£ T -l { T a <+oo}) 
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from which we extract (3.5). The case a < — 1 can be treated exactly in the same way. Finally, 
(3.6) comes from (3.4) by choosing £ = 1 therein and remarking that 1 — Apq = (p — q) 2 . 



Notice the following equality which will be further used: 



P{T_ Q = n} = (J j ¥{T a = n}. (3.7) 
Corollary 10 The distribution of To is given as follows: for even n> 1, 

^->-(;)^ 

and this probability vanishes for odd values of n. Moreover, the probability of eventually returning 
to the origin is given by F{T < +00} = 1 — \p — q\. 

PROOF. By using the Markov property of the random walk as well as the translation invariance, 
we easily see that 

P{T = n} = P{5i = 1} P{T_! = n - 1} + P{Si = -1} P{Ti = n - 1}. 
From (3.7), we observe that pP{T_i = n - 1} = <?P{Ti = n - 1} and then 

, { r = „, = 29 P { r I = „- 1) ^("- 1 )Ml = (;)Ml. 

Thanks to the previous calculations and (3.4) applied to a = 1, we can write out the generating 
function of T : 

00 00 
E (e°l {To<+oc} ) = ]TP{T = n}C = 2 9 £>{T 1 = n - 1}C 



For £ = 1, we immediately find the value of P{T < +00}. 

Remark 11 In the case of the symmetric random walk (p = q = \), the results of this subsection 
take the remarkable following form: 

P{T " = = S (i) if ° € Z V { ° } ' P<T ° = " } = (n-'l)2» (|) ' P{T ° < +~) = 1 ° e 
3.2. Probabilistic interpretation of (1.7) and (1.8) 



By rewriting the distributions of S n and T a (formulas (3.1) and (3.5)) in terms of the Gamma 
function we have that 



Let /x and ^ be positive integers. From (3.8), with i = v — 1 and n = 2r — 2k + v — 1, we obtain 

{r-k)\T{r-k + u) - P q nS2r-2k +v -i-v-l}, 

and, for i = li + v — 1 and n = 2r + ix + v— 1, 

^V^A = P-^^ P(Wi = // + 1> — l}. 
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From (3.9), with a = /x and j — 2k + li, we also have 

F(2fc + M ) _ y-^ q - k _ 2 

Therefore, formula (1.7) can be rewritten as a Darling-Siegert-type relationship (which extends 
formula (4.7) of De Gregorio and Orsinghcr (2007)): 

r 

= 2k+Li} F{S 2r -2k+,-i = v - 1} = P{S 2r+iU+i ,_i = M + ^ - !}• 

k=0 

Similarly, by replacing successively in (3.9) (j, a) by (2k + /x, li), (2r — 2k+v 1 v) and (2r + n + is, n + is) 
we obtain the following expressions 

Formula (1.8) now reads 

r 

^2 P{T M = 2k + /j,} ¥{T U = 2r-2k + v} = P{T M+ „ = 2r + ^ + v}, 

k=0 

which provides a probabilistic interpretation of (1.8). 



3.3. Conditional aspect of (1.7) and (l.l 



Let us divide equation (1.7) by its right-hand side. This yields 



,A r! r(2fc + li) T{2r - 2k + v) T(r + /x + v) 

^ kl (r-k)\ T(2r + Li + u) r(fc + ~jx + 1) T(r - k + v) M ~ ' 



k=0 



that is to say, by means of the binomial coefficients and changing v into v+ 1, that for any integers 
ix > 1 and v > 0, 

r (r\(r+/j,+v\ 

E \k)\ fc+p I M _ 

k=0 \ 2k+n ) ZK + V 



Then, we can observe that the family of numbers {p k v ' T )o<k<r defined by p k ,v ' r — /LV^v 2 k+u 
is a probability distribution. Thanks to (3.1) and (3.5), we recognize that 



„ , , P{T M = 2k + fi} P{S 2r - 2fc+t , = v} 



P{5 , 2r+ M +„ = A» + V } 

Next, by the Markov property and the translation invariance, we have 

P{T M = 2k + fi} P{S 2r -2k+„ = v} = = 2k + fi} ¥{S 2r+ ^ +v - S 2k+f , = v} 

= P{T p = 2k + fi, S 2r+lx+1/ - S 2k +^ = y} = P{T M = 2k + ft, S 2r+Al+ „ = fi + v} 

which leads to the following probabilistic interpretation for p k ' u,r - 

pf = P i T M = 2fc + V\S2r+»+v = M + V}- 

This means that {p k ,u,r )o<k< r is the distribution of the first hitting time of level fi for the random 
walk pinned at the extremity fi + v at time 2r + fi + v. Of course, this interpretation can also be 
directly derived from (1.8). 

Let us consider the particular case fi = 1 and v = and set q r k = p k _'[ 1 for 1 < k < r. We have 

q r k = P{7\ = 2k- l\S 2r -i = 1}. (3.10) 

We also have 

(? £=P{T_ 1 =2fc-l|S 2r _ 1 = -l}. (3.11) 

Indeed, by (3.2) and (3.7), 

|P{T 1 ^2fc-l}P{S 2r _ 2fc = 0} 
~ fP{5 2r _ 1 = l} 

= P{T 1 =2fc-l|S 2r _ 1 = l} = (? £. 

Actually, our aim is to check that 

ql = P{T = 2k\S 2r = 0}. (3.12) 

In this form, we can see that the probability distribution {q k )i<k< r is nothing but that of the first 

return time to for the bridge of length 2r of the random walk. This distribution appeared in De 

Gregorio and Orsingher (2007) in the symmetric case p — q = | and was a consequence of (1.8) 

whereas, here, we use (1.7) instead. Let us now verify (3.12). For this, we write 

P{T = 2k, S 2r = 0} = P{T = 2k} F{S 2r - 2k = 0} 

= (p¥{T-i = 2fc - 1} + q¥{Ti = 2ft - l})P{S 2r - 2k = 0} 
= pP{r_! = 2ft - 1, S* 2r -i = -1} + qPiTi = 2ft - 1, S 2r _i = 1} 
= pP{T_! =2fc-l|5 2r _i = -l}P{5 2r _i = -1} 
+ qP{T 1 = 2k- l\S 2r -i = l}P{S 2r _i = 1}. 



By using (3.10), (3.11), we have P{T = 2k,S 2r = 0} = q r k ( P P{S 2r -i = -1} + ?P{& r -i = 1}) = 
q r k Y{S2r — 0}, which immediately proves (3.12). 



3.4. Probabilistic interpretation of (1.9) 

In the spirit of this paper, we give a probabilistic interpretation for the simple formula (1.9). 
For this, in view of (2.2), we have the equality «! r|/o = 1 where p e (0, 1) is fixed and 

q = 1 — p. This suggests to introduce a family of random variables (T^)^ eN \^, with probability 
distribution 

This is the well-known negative binomial distribution with parameters (/J,,p). Formula (1.9) implies 
that this family satisfies, for any r € N, the relationship 

r 

P{T M = k} ¥{T V =r-k} = P{T^ + „ = r}. 

fc=0 

This formula can be rephrased by introducing an independent copy (T^) M>0 of the family (T M ) AI >o 
as P{T M +T V = r} = P{T AJ+1/ = r}, for any r G N. In words, for any ^, v > 0, the random variables 
T ^ + T v and T have the same probability distribution. If we choose the random variables 
(?V)ai>o independent, then, for any different /i, v > 0, T ^ + T v and X ' M+ „ have the same laws. 

Another interpretation consists of relating the family (T M ) M >o to the family of passage times for 
a certain random walk (S' n ) n &>i, which we are now going to describe. We can rewrite (1.9) in the 
following form: for positive integers \i, v, 

This formula suggests to introduce the following family of random variables. Fix p £ (0,1), q = 
1—p and let (^) MeN \{ j be the Pascal random variable with parameters (fi,p) (this is the negative 
binomial distribution translated by y): 

P{T; =n}={^~_ J^g""", n > (i. (3.14) 

Formula (3.13) reads 

r 

£>{T; = fc + ^{Tl = r - k + v) = P{T; + , y = r + l i + y}. 

Our goal is to find a certain process (<S^)„ eN living on N that obeys the strong Markov property, 
that is translation invariant, and such that for any positive integer /i we have = min{n > 1 : 
S' = /j,} (in order to have S' T , = n) . In view of Proposition 7, for any a € N, the generating function 
G'(£,a) defined, for |£| < 1, by G'(£,a) = Y,n=o F i S n = a R" satisfies the following relation 

G'(Z,a)=G'(Z,0)E(e°). (3.15) 



Above, the generating function of T' a is explicitly given, with the aid of (2.2), by 

e(^.) = £ (; : J)^— r = £ ( n +! ; ^vr^ = (^y " (3.16) 



and we have 



CK,0)~ (3-17) 
1 - 

since, by observing that we obviously must have J2T=o ^{^'n = a} = 1 and by referring again to 
(2-2), 

^ oo oo oo oo oo / j. \ a 1 



,_, E^ = EE p K=^ B = E G '(^) = G ^ )E t^t =guo) , 

n =0 n=0o=0 a=0 a=0 \ y ^ / 

By plugging (3.16) and (3.17) into (3.15) and using (2.2), we obtain that 

g'(^) = (T^fpr = wrE ( n+ a a yc - W"EQrT-» = E Qpv-r 

which implies that 

P{S' n =a} = (^jp a q n - a , 0<a<n. (3.18) 

So, we recognize the binomial distribution of parameters (n,p). We deduce from this discussion 
that the process we were seeking is the non-decreasing random walk (S' n ) n eN where 5^ = J2k=o ^'k 
with Xq = and such that the X' k , k > 1, are independent random variables taking the value 1 or 
with probability p or q, respectively. 

Conversely, we can easily check that this random walk satisfies the three assumptions (strong 
Markov property, translation invariance, discrete continuity) made in the above analysis. 

This random walk plays an important role in game theory. Indeed, if we toss n times a coin and 
call "head" a success, S' n is the number of successes among the n tosses while T' a is the minimal 
number of tosses for obtaining exactly a successes. 

We have interpreted (1.9) by means only of the variables (7^) aeN \{ }. Since we have identified 
the underlying random walk (S' n ) n eN associated with (T„) o eN\{0}> ^1 us focus ourselves on the 
Darling-Siegert relationship connecting both families. We have F{S' n = a + b} = Xwfc=a^{^o = 
k}^{S' n _ k = b} which, by (3.14) and (3.18) supplies the simple combinatorial identity 

(k — 1\ fn — k\ ( n 

Rewriting this latter, by setting n = r + a + b, as 



b J \a + b 



^/fc + a-l\/r-fc + &\ _ fr + a + b^ 
h\ a ~ 1 b )-{ a + b 

we retrieve, by translating one of the parameters a and b, formula (3.13). We then conclude that 



both Darling-Siegert-type relationships related to (S' n ) ne ^ together with (T^) agN \{ } 



i-6 



F{s' n = a +b} = j2nK = k} nsu = b}, nK +b = »} = E P ^ T - = fc > w = n - fc > 



k— a k—a 



lead, in this case, to the same combinatorial identity, contrarily to the foregoing case. This fact 
is in good accordance with the observation that both identities (1.5) and (1.6) lead to the same 
combinatorial identity (1.9). 

Remark 12 The particular relationship P{T a = n} = ^ P{S n = a}, for any a e N \ {0} 7 which 
holds for both random walks considered in this paper, shows that one Darling- Siegert-type formula 
implies the other one. Indeed, starting from Efc=o ^{^o = k}P{S n -k = b} = P{S n = a + b}, 
we obtain ELo( n _ k )H T a = k}P{T b = n - k} = ¥{T a+b = n}. By making the change of 
index k — > n — k and by interchanging the roles of a and b, this can be written as Efc=o ^¥{T a = 
k}P{Tb = n — k} = ^pjP{T a+ b = n}. Now, adding the two foregoing equalities leads to formula 
ELo HTa = k} P{T b = n-k} = V{T a+b = n}. 
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